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ABSTRACT 



Current theories of galaxy formation have tended to focus on hierarchical structure 

formation, which is the most likely scenario for cosmological models with lots of 

power at small scales (e.g. standard cold dark matter). Models with little small scale 

power lead to scenarios closer to spherical collapse. Recently favored power spectra 

■ (e.g. CDM+A) lie somewhere in between suggesting that both types of processes are 

CN ■ important and may vary over time due to gaseous reheating. From this viewpoint this 

paper explores a very simple inside out scenario for galaxy formation. This scenario is 

a natural result of synthesizing earlier work on DM halos, spherical collapse, and gas 

redistribution via angular momentum. Although, this model is highly simplified and 

is not designed to accurately describe the detailed formation of any individual galaxy, 

it does (by design) predict the overall features of galaxies. In addition, old bulges 

and young disks are an almost unavoidable result of these very simple models. This 

> ■ scenario may provide a useful framework for both observers and theoreticians to think 

c/3 ■ 

& . about galaxy formation. 

> 

X 

1. Introduction 

Galaxy formation theories can be caricatured by two models: hierarchical clustering and 
spherical collapse; with reality most likely lying somewhere in between. The degree of hierarchical 
clustering is strongly dependent upon upon the amount of power on small scales. Cosmological 
models with lots of small scale power (e.g., standard CDM) will be more hierarchical, while many 
of the cosmological models currently under consideration have power spectra with less small scale 
power (e.g., mixed dark matter, tilted CDM+A) ( |Borgani et al 1997 ) and will be less hierarchical. 



The thermal history of the gas can also play an important role. At early epochs the gas 
temperature will be low, which allows small lumps to form and be accreted hierarchically. At later 
epochs, the IGM heats up to ~ 10 4 °K due to the first generation of stars and Quasar emission. 
Gas at this temperature will be driven out of halos with circular velocities less than 30 km s _1 



( Kepner, Babul &: Spergel 1997 ). Accretion of these small lumps at z ~ 3 may not take place 



hierarchically but will tend towards spherical infall as the puffed up lumps of gas encounter gas in 
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the halo that has been heated by shocks and radiation from a variety of sources: stars, supernova, 
and possibly a super massive black hole. In reality, the gas in galaxy formation is most likely in 
multiple phases (like the present day galaxy). 

Exploring both hierarchical and spherical collapse aspects of galaxy formation is best 
accomplished with a variety of techniques. N-body simulations are best for exploring hierarchical 
structure formation. Detailed simulations of the formation of individual galaxies can now be 
performed with an ever increasing list of physical processes. Some of the important results of 
simulations are the approximate correctness of Press-Schechter theory for describing merging 



histories ( Lacey &; Cole 1994 ), that dark matter (DM) halos have a similar form over a wide range 
of scales ( Navarro, Frenk fc White 1996|) , and the nature of the Lya forest ( Miralda-Escude et al 
1996, Ma et al 1997, Zhang et al 1997] ). In addition, simulations that include gas dynamics show 



development of disks as a result of tidal torques ( [Navarro Sz Steinmetz 1997 , Steinmetz &; Ewald 



1995) 



Analytic methods lend themselves to the understanding of spherical collapse. The original 
analytic work on spherical collapse (Gunn & Gott 1972, Fillmore Sz Goldrcich 1984, Bertschingei 



1985) laid the groundwork for understanding virialization, self-similar collapse, the establishment 



of near isothermal dark matter halos, and secondary infall. More recent analytic work has been 



particularly useful in providing intuition on the formation of disk galaxies ( palcanton, Spergel 
fe Summers 19971 , pVlo, Mao fe White 1997|) , an area that is particularly challenging for N-body 
simulations due to the large dynamic range requirements. These recent analytic papers, built upon 
the earlier work of Mestel 1963| , Fall fe Efstathiou 1980| , and Gunn 198% , indicates that allowing 
gas with typical angular momentum distributions to settle into typical DM halos produces disk 
galaxies that reproduce a wide variety of observations from faint surface brightness objects to 



damped Lya systems. The fundamental assumptions of Dalcanton, Spergel & Summers 1997 
and |Mo, Mao fc White 1997 , which will also be adopted in this paper, are that (1) there is little 
angular momentum transport, and (2) the angular momentum eventually halts the collapse of the 
gas, resulting in a roughly flat rotation curve. 

Spherically symmetric analytic models have a limited ability to address highly non-linear 
phenomena such as hydrodynamic shocks. Recent spherical numerical simulations of infall into 
galaxies indicate that gas can undergo a shock before settling into the DM potential ( |Thoul fc 



Weinberg 1995| ). Such shocks are potentially an important process for slowing down and heating 
infalling gas during galaxy formation. 

The goal of this paper is to add to the analytical understanding of galaxy formation by 
connecting the earlier spherical collapse models with the more recent work on redistribution by 
angular momentum and infall shocks. The approach is quite simple. Given a known approximate 
initial state — some kind of primordial "bump" ( |Bardeen et al 1986 ), and a known approximate 
final state of the DM halo and the gaseous disk, completing the picture requires determining the 
time, mass and length scales for the transformation of the DM and the gas from the initial state to 
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the final state. The collapse and virialization of spherically symmetric DM halos gives an estimate 
of the scales involved in creating a DM halo. The scales for gaseous collapse can likewise be set by 



the evolution of an outwardly moving shock (Thoul & Weinberg 1995). 



Synthesizing these two components results in a single physically motivated scenario which, by 
construction, reproduces many of the observable quantities of the final states of galaxies, but also 
makes specific predications about gas evolution during galaxy formation. Such a scenario can be 
a useful perspective for both observers and theoreticians by clarifying the important quantities 
to be measured and the additional physical processes that need to be included to explain various 
phenomena. 

In §2 a sketch of the inside-out galaxy formation scenario is presented along with an overview 
of the underlying physics. In §3 the physical model is presented in detail. §4 discusses the results 
of these calculations. §5 gives conclusions and plans for further work. 



2. Inside-out Galaxy Formation 

Any exploration of galaxy formation begins with a list of ingredients. The two main 
components are dark matter and gas. In this paper it is assumed that the dark matter is cold and 
collisionless. The gas can behave in a variety of ways depending upon its thermal history. If the 
gas is dense it can radiate away heat acquired in shocks and will tend to behave in a relatively 
cold manner. If the gas is diffuse it may be heated by external radiation or shocks. Recent 
simulations of infalling gas with radiative cooling indicates that such shocks occur in a variety of 
galactic potentials ( [Thoul &: Weinberg 1995| ). Furthermore, these simulations suggest that the 
initial radius at which a mass shell is shocked is similar in both the adiabatic and non-adiabatic 
cases. Thus the simplest approach is to use a adiabatic numerical approach to calculate the initial 
shock radius and then use analytic non-adiabatic methods to estimate the post-shock behavior. In 
reality, the gas in galaxy formation is most likely in multiple phases (like the present day galaxy). 
Of course, the limits of the adiabatic approximation for computing the shock radius must be kept 
in mind in interpreting any subsequent results. 

The inside-out scenario that results from combining spherical collapse and angular momentum 
redistribution proceeds as follows. Initially the gas and dark matter are coupled and expand 
with the Hubble flow. The innermost shells turn around first. The dark matter virializes, and 
the gas is shock heated to the virial temperature. The cooling time for these inner shells is very 
short and some of the gas may quickly condense into cold lumps and form bulge stars. Later, the 
outermost shells turn around and their gas is shock heated. However, for these outer shells, the 
higher virial temperatures and lower gas densities result in longer cooling times; so that as the gas 
cools and falls to its corresponding angular momentum radius it is shocked again and forms stars. 
Higher angular momentum material ends up outside the bulge and forms the disk. Lower angular 
momentum material falls into the bulge, which continues to grow slowly, but may be halted by 
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energy input from supernova and/or a massive black hole that can initiate a wind that blows out 
the gas. 

The two components of the above scenario are DM and gas. In the next two sub-sections the 
specific evolution of each component is elaborated. In addition, some of the concepts that form 
the basis of the physical model will also be introduced. 

2.1. DM Evolution 

Typically, a spherical dark matter halo is quantified by a characteristic mass and radius 
(-^haio> r haio)) which can be translated into a circular velocity and virialization redshift (Vcn-c^riai) 
by assuming a value for the overdensity at z v i r i a i (e.g., <5 v i r i a i ~ 179). Simple estimates of the 
number of halos of a given mass at each redshift can be made via the Press-Schecter formalism. 
These estimates provide useful insights into the size distribution of DM halos and are in 
approximate agreement with galaxy catalogs and N-body simulations. 

As primordial DM halos initially expand and contract they are torqued by neighboring 
density peaks leaving some galaxies with significant amounts of angular momentum. The amount 
of angular momentum is characterized by the spin parameter A sp i n . Galaxies with larger values 
°f ^spin form disks. Thus three parameters would seem to describe the gross characteristics of 
galaxies (Mh a io,nialo,A S pin) or equivalently (V c i rc 

For an idealized spherically symmetric DM halo, the evolution is qualitatively quite simple. 
Each mass shell in a primordial "bump" proceeds through Hubble expansion, turn around, 
virialization and relaxation into a final quasi-stationary DM profile. This process is characterized 
by the collapse time i C oiiapse and the maximum radius r max of each bound DM shell. 

The initial bump can be estimated via various theoretical arguments and the result can 
usually be fit by an approximately Gaussian profile. The final DM profiles have been extensively 
studied both observationally and via N-body simulations and a variety of universal profiles have 
been proposed. These differ in detail but produce similar rotation curves, V c i rc (r), over the 
observationally accessible radii (see Figure |]) . The process of virialization has also been extensively 
studied in simulations. There is no simple way to characterize this process. Fortunately, for 
spherical collapse this process seems to take place fairly quickly and in a manner which tends to 
erase the details of the initial distribution ( [Hoffman 1988| ). 

An asymptotically correct model of the DM evolution can be constructed by marrying the 
initial expansion and turn around of a pre-determined initial profile with a pre-determined final 
profile by simply stopping each infalling mass shell at the radius set by the final profile. Such a 
model — while completely sidestepping the important role of virialization and hierarchical structure 
formation — provides a simple framework for describing the DM potential and its effect on the gas. 
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2.2. Gaseous Evolution 

For the most part, each gas shell of the galaxy evolves passively under the gravitational 
influence of the dark matter until it collapses. In the case of adiabatic gas, the inner shells 
will be compressed until there is sufficient pressure to stop the infall and create a shock. This 
process is characterized by the shock time t s hock an d shock radius r snoc k of each mass shell which 
define the pertinent time and length scales of the gas in the halo. Inside the shock, the gas will 
be approximately in isothermal hydrostatic equilibrium with a temperature given by the DM 
potential at r snoc k. Outside the shock, the gas will be cooler and significantly less dense. In a real 
multi-phase gas, r s h oc k can be thought of as an approximate upper limit for the radius at which 
cooler clumps of gas begin to feel the effects of the external environment. The "true" radius will 
be somewhat smaller to account for the fraction of gas that has condensed into cold lumps. Some 
simulations indicate that early on this can be a large fraction of the gas ( |Navarro fc Steinmetz 



1997 ). Later on the ratio of cold gas to shocked gas can be as low as 20% ( [Thoul fc Weinberg 



1995|) , suggesting that the adiabatic value of r sn0 ck becomes more accurate at later times. 



The main effects of the shock are twofold. First, the dissipation of a significant amount of 
inward velocity results in the slow spiraling of any cool clumps towards their corresponding angular 
momentum radii. Second, a large increase in density leads to a dramatic drop in the cooling time 
*cooi oc /Jgag. If the cooling time is short compared to the dynamical time t^ yn = ^Kirc' then even 
the clumps in the hot gas will condense and may undergo rapid star formation, otherwise they will 
spiral inwards until they hit the disk. The initial gas mass and radius of this region, which may 
correspond to an initial bulge, can be computed by setting idyn = *cooi- Subsequent mass shells 
will split their material between the bulge region and the disk depending upon the distribution 
of specific angular momentum within each mass shell. The bulge will continue to grow slowly. 
However, this process can be halted if there is sufficient energy input from supernovae and/or a 
massive black hole to drive a wind. The disk will continue to grow outward as long as there is 
infalling material to supply it. 

Asymptotically redistributing the gas according to its angular momentum produces disks in 



agreement with a wide variety of observations ( Dalcanton, Spergel fc Summers 1997 , Mo, Mao & 



White 1997] ). Thus, as with the DM halo, there is a fairly good picture of the early and late stages 



of the gas. The evolving adiabatic shock is a means for setting the time, length and mass scales of 
the evolution of the gas between these states. 



3. Physical Model 

The previous section gave a broad overview of an inside-out scenario for galaxy formation 
and hinted at the underlying physical model used to motivate it. In this section the details are 
presented. First, the final DM halo is described. Second, the evolution of the DM halo from a 
primordial bump is discussed. Third, the evolution of the shock radius of the gas is presented. 
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Fourth, the initial formation of the bulge is computed. Fifth, the infall of material onto the bulge 
and the disk is derived. Sixth, the processes that may lead to a wind are discussed. 



3.1. Final DM halo 

Let the final dark matter halo be specified by two parameters: the halo radius and halo mass, 
which can be translated into the the circular velocity V c \ rc and the virialization redshift z V mai if it 
is assumed that the overdensity at virialization is <5 v iriai ~ 179 (Gunn & Gott 1972) 



T/2 / \ GM ha i 47T 3 

^circl r haloj — , ^ r h a lo<AdrialPcrit l^virial ) ~ -Mhalo> UJ 

^halo o 

where the mean density is given by the usual expressions for a f2 = 1 CDM cosmological model: 
Pcrit(z) = (1 + zfp® rit , 67rGpcrit*;Liverse = ^ Universe = 2 /3ff . The distribution within the 
DM halo can be any reasonable profile (e.g. Hernquist 199C| , [Burkert 1995 , |Navarro, Frenk fc 



White 1996, see Figure hi), which satisfies the above equations at r^aio- The simplest choice is an 



isothermal sphere for which: 



M DM (r) = Mhal ° r , V circ (r) = V circ (r ha i ) , pdm(^) = Mhal ° _ i_ (2) 
^haio ^haio 4vrr^ 



3.2. Primordial bump 

The simplest description of the state of a DM halo at some early epoch (corresponding to 
redshift Zinit) is a spherical "bump" in the density field /0dm(?"> ^mit)- A typical initial density 
profile might be (see Figure ^): 

PDM(r, Zinit) = Pcrit(zinit) + p e~ r2 /r ^° {zinit) (3) 



The dynamics of each mass shell in the halo are governed by the interior overdensity: 

M DM (r 

^-r 3 p C nt(Zmit) 

where Mdm(?") = Jq ^r' 2 puM(r')dr' . 5dm determines the Hubble flow for each point at Zmit via 
the linear theory estimate (Thoul & Weinberg 1995) 

VDM(r, Zinit) = rH(z init )[l - -<5 DM (r)]. (5) 
The above estimate will be accurate while 5dm « 1. 

Having specified the initial position, velocity and time of each mass shell the subsequent 
evolution is governed by 

o t a ■ a 2r i a 7Ts[n6 rmax fa\ 
2ir- = 0-sin0, = l-cos0, i>dm = — — , (6) 

'•collapse '"max J- COS (7 I co n a p SC 
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where 

1 __ 1 v ( r ) i \ _ / r max / 7 \ 

W0 = r"2GM DM (r) ' tcollapse(r) = ^M(0 (?j 

The above description of the trajectories holds as long as there are no shell crossings. 
Unfortunately, there is no elegant procedure to ensure collapsing shells virialize into a desired 
quasi-stationary DM halo. The simplest method, which is adopted here, is to stop each mass 
shell when it reaches a pre-determined final radius. It remains to be seen whether or not this 
approach provides a sufficiently accurate picture of the DM potential from the point of view of the 
gas. Figure ^ shows the trajectories of various mass shells when this method is employed with an 
isothermal profile. 



3.3. Shock Radius 



Detailed spherical collapse calculations ( [Thou! fc Weinberg 1995| ; |Thoul fc Weinberg 1996|) 



indicate that gas falling into a variety of galaxy potentials undergoes a shock. These simulations 
were conducted for both adiabatic and non-adiabatic gas (using standard cooling functions) in a 
fully dynamical DM potential. One important observation of these simulations is that the initial 
shock radius in both adiabatic and non-adiabatic gas is similar. Howver, after the initial shock 
the gas cannot be treated adiabatically. While a non-adiabatic model in a fully dynamic DM 
potential provides a richer description of galaxy formation it is more difficult to develop a simple 
relations for the shock radius from such simulations. This would suggest that the simplest means 
for obtaining an expression of the initial shock radius is to examine adiabatic gas falling into a 
deterministically evolving DM potential. After obtaining this relation, a non-adiabatic analytic 
approach is adopted in the next sub-section to look at the post-shock behaviour of the gas. 

The adiabatic evolution of the gas within the above evolving DM potential is governed by the 
Lagrangian equations: 

-^-M gas = 47rr 2 /9 gas , 
± v 2 dP gas GM DM (r) 



dt gas " dM gas r2 

d ^ -fgas d 

Jt u ^ - 7^dt p ^' 

Pgas = (7 - l)Pgas% as i (8) 

where the self-gravity of the gas has been ignored because it becomes important only in the 
inner regions of the bulge and the disk, well inside the shock. For the relavent initial conditions 
(Pgas(f) = ^baryonPDM(^)) v sas (r) = vdm the- global behavior of the solutions to the above 
equations are quite similar: the gas expands and contracts with the DM until the pressure gradient 
stops the infall and causes a shock to form. Inside the shock the gas is hot and dense, while outside 



- 8 - 



the shock it is cool and less dense. Thus, the evolution of the shock front describes, to first order, 
the overall behavior of the gas. The shock occurs when the infall speed of a shell approximately 
equals the thermal speed of the gas interior to it. If the gas is assumed to follow the DM until this 
point, and that inside the shock the gas is at a temperature corresponding to the circular velocity 
of the DM halo, then the shock will occur when: 

^circ ('"shock) ^gas (A^shock) ~ ^DMv^shock/^baryoii) ^shock) > (9) 

where M shock = M gas (r shock ) = 0, ha _ Tyon M DM (r s hock)- Figure g shows plots of 

'"shock VS. M shock for 

the four halos shown in Figure [l] obtained from solving the above equation and compares it to 
that obtained using a spherically symmetric Lagrangian hydro-dynamics code similar to that of 
Thoul fc Weinberg 1995| . Both results obey a simple power law 

/^shock 

, fthock « 0.82. (10) 



'"shock / ^shock 



'"halo \ ^baryon-^halo / 

In the next two sections this result is used to estimate the gas properties of the bulge and the disk. 



3.4. Initial Bulge 

Inside the shock the gas is approximately isothermal with a temperature given by 
Tshock = arKh-c^shock)- In addition, the gas will be in hydrostatic equilibrium. In the case of a 
DM halo with a flat rotation curve, the density of the gas just inside the shock is 

Pshock — Pgas ('"shock) ~ 4^ 3 ; (1-Q 

3— ctr ^shock 

If the cooling time is short compared to the dynamical time, then clumps in the shocked gas will 
condense and form the first bulge stars, otherwise they will spiral inwards until they hit the disk. 
The boundary between gas that cools quickly and gas that hits the disk can be computed by 
setting t dyn (r shock ) = t coo \(r 

shock ) • 

ft r shock _ Ashock^B^shock ^ (-10} 

~ (7 - l)/i A(T shock ) Ps 2 hock ' 1 j 

where A(T) is the primordial optically thin cooling function ( [Thoul &: Weinberg 1995| ), [i is the 
mean molecular weight, ks is Boltzman's constant and 7 = f . Manipulating the above expression 
gives 

^4hock _ 8 pfceTsho^Vrirc . . 

dock " 3 - or (7 - l)A(Tshock) 1 j 

All the terms on the right can be obtained directly from the parameters specifying the DM halo. 
If M s hock an d '"shock are related by a simple power law, then the initial mass and radius of the 
bulge can be estimated. Figure I shows rj^ and M^ gc for four DM halos MS M function of V c i TC 
for cut = 2. 
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3.5. Disk Formation 

Subsequent mass shells will split their material between the bulge and the disk depending 
upon the distribution of specific angular momentum within each mass shell. The simplest model 
for the initial angular momentum distribution is to assume that each shell is in solid body rotation 
with frequency ujq when r = r max , where ujq = A sp i n V^i rc /rh a i . This model is consistent with the 
idea that a typical proto-galaxy acquires the majority of its angular momentum from an single 
encounter with a similar sized object when it is at maximum expansion. In the bump initial 
conditions used here the inner shells will turn around before the outer shells. If the external 
torquing field changes in time, then the shells will have different angular momentum histories. 
This is most likely the case in the real galaxies. However, because solid body rotation does an 
adequate job of explaining the current state of disk galaxies ( palcanton, Spergel fc Summers 1997 , 



Mo, Mao & White 1997), it is convenient to maintain this simplifying assumption. 



The distribution of specific angular momentum within a rotating shell is not the same. Gas 
at the equator will have more angular momentum than gas at the poles. This distribution of 
specific angular momentum within a gas shell at a time when the shell is at a radius r gas is 
igas = (-Rga S /Vgas) 2 .?max, where j max = w r^ ax and R gas < r gas is the projected radius. The specific 
angular momentum of each radius in the DM halo is j'dm = V^; rc (i?DM)-RDM- The final distribution 
of the gas can be computed by setting j gas = j'dm 

M gas ( J R DM ) = jl - ^I-Jdm^dm)/^! M*°*, (14) 

The majority of the mass of each shell ends up near the maximum radius Jdm(-Rdm) = Jmaxi which 
for an isothermal halo is Rbm = Jmax/K:irc- The final distribution of the disk can be computed by 
summing the contributions of each gas shell. Detailed calculations of this type reproduce a wide 



range of observed properties (Palcanton, Spergel fc Summers 1997, Mo, Mao fc White 1997). 



The gas takes approximately idyn( r shock) to settle into the above distribution after it has been 
shocked. Because the majority of each gas shell falls near the maximum angular momentum radius, 
the disk grows outward according to r^(t shoc k(Af s hock) + t dyll (M shock )) = j ma , x (M shock )/V circ . 
This outer edge will grow as long as there is infalling material. 

In contrast, the bulge continues to grow by slowly accumulating the low specific angular 
momentum material of many gas shells. The total amount of material falling into the bulge can 
computed by integrating the shocked shells 

-Wbulge(ishock(M shock ) + tdyn(M shock )) = / I 1 - yjl - jDM^bulgc) / jmax(M s ' hock ) > dM s ' hock 

(15) 

In general, t dy n + ishock ~ Collapse, so for any halo: 



M bulgc (t colla pse(Af shock )) = / MSh ° Ck ( 1 - Jl - - } dM & ' hock (16) 

JO { V V / circ(^halojAspinrmax(JW s hock) J 
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The gas in the disk is simply M<ji s k = Mohock — Afb u i ge - Figure |6| shows the evolution of the infall 
rates Mbuige and Mdisk as a function of redshift for four different DM halos. Interestingly, this very 
simple picture produces old bulges and young disks quite naturally, with reasonable values for the 
total mass and infall rates. 



3.6. Blowout 

The bulge will continue to grow slowly without additional sources of energy to heat the 
gas and drive a wind. There are three potential sources of energy input: Supernovae Type I, 
Supernovae Type II and a massive black hole. A realistic treatment of the feedback due to any 
of these processes is a challenging undertaking and is fundamentally limited by an incomplete 
understanding of several processes (e.g., star formation, supernova IMF, black hole formation, 
etc ...). Here only the most simple estimates are attempted. Fundamentally, to drive a wind the 
energy input into the bulge gas must overcome the escape velocity of the infalling gas: 

-^bulgeKsc < max { Q; SNI-t'SNI, CKSNII-^SNII, OSNI-^Bh} (17) 

where V esc is the typical escape velocity of gas in the bulge, and q;sni,snii,bh is the fraction of the 
total energy input of each type which is usefully deposited into the gas in such a way as to drive 
a wind, i.e. a parameterization of our ignorance about converting luminosities into winds. The 
energy inputs can be estimated as follows 

-^SNI = eSNI-^bulge^SNlAsNI (18) 
-kSNII = esNIlMbulgeVgNn (19) 

Lbh = ebh-Mbhc 2 ~ L cdd = M B uc 2 /t edd (20) 

where , t edd = 47r CT J^ H ~ 5 x 10 9 yr. 

The equation for Lsni assumes the progenitors of Type I SN (really Type la) are white 
dwarfs or some other compact stellar remnant which should be proportional to the total bulge 



population. If the average rate of SNI is one per hundred years per 10 1U M ( pain et al 1996|) and 



each SNI ejects approximately 1 Mq of material, then setting esNi ~ 0.01, Vsni ~ 10 4 km s 1 and 
tsNi ~ 10 10 yr will give approximately the correct rates and energies for SNI. 

The equation for Lsnii assumes the progenitors of Type II SN (really Type Ib,c and Type II) 
are massive short lived stars whose rate should be proportional to the star formation rate, which 
in turn is roughly proportional to the infall rate. If the average rate of SNII is once per hundred 
years ( Cappellaro et al 1997] ) in galaxies with star formation rates and infall rates of on the order 
of one solar mass per year, and each SNII also ejects ~ 1 Mq of material, then setting esNii ~ 0.01 
and Vsni ~ 10 4 km s _1 will give approximately the correct rates and energies for SNII. 



The equation for Lbh assumes the massive black hole radiates a fraction of the gravitational 
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energy of the material falling onto it. Furthermore, it assumes that the luminosity and subsequently 
the infall rate is constrained by the Eddington luminosity of the black hole. 

Using these definitions it is possible to try and guess under what conditions each type of 
energy input will drive a wind. The blowout condition for Type II Supernova is the simplest 

i i 

V esc < ^SNII«SNII e SNII ' 

(21) 

The value of cksnii is difficult to estimate, but can be obtained by setting V esc ~ 100 km s _1 the 
equivalent critical value arrived at in more detailed calculations flDekel fc Silk 198(f ), in which case 
cksnii ~ 0.01. 

For Type I Supernova the blowout condition is 

-^bulgeKsc < a SNI e SNlMmlge^SNlAsNI > (22) 

which can be rewritten in terms of the bulge growth time Mb u i ge /Mbuige 



/ "SNieSNI-Mbulge 



V csc < v sm ■ ( 23 ) 

\ tSNPWbulgc / 

The above condition is plotted in Figure |^ (along with the SNII condition) assuming the 
QfSNi = cksnii- 

This oversimplified model does not readily produce a black hole blowout condition, but it is 
possible to estimate the ratio of the black hole mass to the bulge. If the black hole accretion is 
limited by the Eddington luminosity, then Mbh = -^bh ex P[*/ e BH^edd]> which sets the time scale 
for blowout tbiowout = fBH^edd ln[MBH/-^BH]- Combining the equation for black hole luminosity 
with the blowout condition gives 

M BH > ^fMbuigc (24) 

If accretion onto the bulge is relatively constant, then Mb u i ge ~ MbuigeAbiowout and the ratio of 
the black hole to the bulge at blowout is 



. y ■ esc/ ~j (25) 

M bulg J,, f «BHeBHln[MBH/M^] ! 



MbH_\ . (Ksc/C) 

X % e J blowout 

which for the characteristic numbers €bh l n [^BH/^BH] ~ 1 an< ^ ^esc ~ 100 km s _1 requires a 
value of aBH > 10 -4 in order to be in agreement with the observed black hole to bulge ratio of 
~ 10~ 3 ( [Faber et al 1996 ). Using this bound for obh we can now estimate the blowout condition 

V„<(*^^) h (**L)* , (26) 
V *edd M hulgc J \M bulge J blQwout 

which is also shown in Figure [7| for obh = 10 -3 . 
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4. Discussion 

The previous sections synthesized an inside-out galaxy formation scenario by combining a 
spherical collapse model with the final redistribution of gas by angular momentum. The aim was 
to interpolate the gaseous evolution from models of the initial and final states of the DM and the 
gas. 

The DM is governed by the specified final halo profile and r max and i C oiiapse> which are set 
by the mass overdensity of the initial profile 5 m . The primary focus here is on the gas, and its 
response to the DM potential. Thus, a highly simplified model is adopted of the post-collapse 
state of each DM halo. Each DM shell stops falling at the radius given by a pre-set final DM 
profile. The validity of this approximation is questionable. Clearly it cannot properly account for 
any hierarchical buildup of material in the halo. If, however, DM halos — in spite of hierarchical 
accretion — maintain a quasi-self-similar density profile as is suggested from N-body simulations 
( Navarro, Frenk White 1996| ), then this asymptotic approach might be OK. 



The state of the gas can be characterized by the evolution of the adiabatic shock in response 
to the above DM potential. For Gaussian initial conditions, r s h oc k and M s i loc k can be related by 
simple power law (Figure ||) for a variety of DM halos. As mentioned earlier the adiabatic value 
of ?"shock(-^shock) is only an upper limit. In a real multi-phase gas, r s h oc k will be smaller due to by 
the fraction of gas residing in cold clumps. This effect will be most pronounced earlier on, but 
during later stages of the infall the adiabatic approximation is in good agreement with similar 
calculations that take cooling into account ( [Thou! fc Weinberg 1995| ). 

Using the power law for the shock evolution, the size of the initial bulge (Figure ||) can 
be estimated by computing the cooling time just inside the shock and setting it equal to the 
dynamical time computed from the DM halo. It is reassuring that the initial bulge mass and 
radius are reasonable and display the correct dependence on the halo parameters, except for the 
Burkert profile, which is only meant for dwarf galaxy halos with V C \ TC < 50 km s _1 . Unfortunately, 
the sensitivity of these estimates to ay and to the details of the DM halo make this approach 
unsuitable for precise calculations. 

After the infalling gas has been shocked, if it does not cool quickly and form stars, then it will 
spiral in to its corresponding angular momentum radius in about a dynamical time. The simplest 
model for the angular momentum of the gas is to assume that when each gas shell reaches its 
maximum radius it is in approximately solid body rotation with a frequency determined by A sp i n . 
What is most interesting is the relative ease with which one gets an old bulge and a young disk 
from this highly idealize picture. The infall rates (Figure |||), which are upper limits, are only a 
factor 10 greater than those seen observationally. The size of the disk vs. the bulge depends upon 
the quantity rbuige^haio^spirn thus there is a large amount of freedom to create diverse objects. 
More complex angular momentum models ( |Ryden fc Gunn 1987 , |Eisenstein &: Loeb 1995] ), which 



take into account the spectrum of initial fluctuations and estimate the statistical distribution of 
angular momentum, could also be used. The overall effect of these models would be to shift and 
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soften the transition between bulge and disk infall. 

The bulge will continue to grow slowly and may be halted if there is sufficient energy input 
from supernova and/or a massive black hole to drive a wind. The disk will continue to grow 
outward as long as their is infalling material to supply it. The simple blowout estimates computed 
in the previous section indicate that all the mechanisms examined have the potential to drive a 
wind which may affect the accretion onto the bulge. 



5. Conclusions and Further Work 

A simple inside-out scenario for galaxy formation is presented, which is a natural result 
of synthesizing earlier work on DM halos, spherical collapse, and redistribution via angular 
momentum. The scenario predicts a straightforward sequence of events beginning with the 
formation of an initial bulge due to the high densities and rapid cooling times that exist when the 
first infalling mass shells are shock heated. After the initial collapse, the bulge continues to grow 
slowly but may be halted by energy input from supernova and/or a massive black hole initiates 
a wind that blows out the gas. A disk forms as the outermost shells collapse. The higher virial 
temperatures and lower gas densities of these shells result in longer cooling times; so that instead 
of quickly forming stars the gas spirals in towards its corresponding angular momentum radius. 

Although, this model is highly simplified and is not designed to accurately describe the 
formation of any individual galaxy, it does (by design) predict the overall features of galaxies. 
In addition, it appears that old bulges and young disks are an almost unavoidable result of this 
very simple model. Furthermore the epochs at which maximum infall into the bulge and inner 



disk take place are roughly consistent with both the high star formation rates flMadau, Pozzetti 



fc Dickinson 199% ) and quasar densities ( Boyle et al 1997 ) observed at z ~ 2. A blowout taking 



place early in the evolution of a galaxy is also consistent with the observation of high velocity 
winds seen at high redshift ( [Pettini et al 1997| ). This scenario also suggests two mechanisms for 
preventing the formation of disk galaxies in the center parts of clusters. First, by encounters which 
disrupt the alignment of the angular momentum vectors of the infalling gas. Second, by the high 
virial temperatures present, which would cause late infall to be absorbed into the ambient cluster 
gas. Of course, these are only two of many possible explanations, but never-the-less they may be 
worth exploring with more detailed simulations. 

These simple calculations suggest an idealized framework in which to attempt the direct 
simulation of disks. Cosmologically based simulations are limited by the need to construct the 
potential from particles, which leads to a significant increase in the complexity and duration of 
the computations. Observing 3D gas collapse into a spherically symmetric time evolving potential 
like the one presented here should be a useful test case of the 3D codes. Further work in this area 
should also focus on exploring the role of the processes that intrinsically cannot be included in 
this simple model: hierarchical structure formation and a multi-phase ISM. 
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Finally, the extremely simple estimates of the blowout condition suggest that Type I and 
Type II Supernovae and super massive black holes all have the potential to initiate an outflow 
and subsequently change the evolution of the bulge. Further calculations along these line using 
more realistic assumptions could prove fruitful in explaining the observed properties of bulges and 
elliptical galaxies. 
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Fig. 1. — The circular velocities for different profiles for a DM halo with V^i rc (rh a i ) = 200 km s _1 , 
-Zvirial = 3, -Mh a io = 10 12 Mq, rhaio = 54 kpc. All the profiles reproduce typical rotation curves over 
the observable range. Hernquist: M oc r 2 /(r + a) 2 , isothermal: Mar, Navarro: M oc log(r 2 + a 2 ), 
Burkert: p oc [(r + a)(r 2 + a 2 )] -1 , where 5a = r^aio- 
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Fig. 2. — Initial density profile and initial mass overdensity profile obtained from a Gaussian bump, 
and a corresponding top-hat profile for the halo given in Figure 1 (zinit = 68, rhaio(^imt) = 16 kpc). 
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Fig. 3. — Trajectories of selected DM shells. The shells are stopped at the radius corresponding 
to the final position for the isothermal halo given in Figure 1. 




Fig. 4. — Shock radius r s hock y s shock mass M s h oc k for the four DM halos given in Figure 1 using 
three methods: (i) Lagrangian hydrodynamics code (ii) infall velocity constraint ^^(rghock) = 
v DM (M shock ) (iii) simple power law (r sho ck/niaio) = (Mshock/^aryonMhaio)^ 5110 *, /3 s hock = 0.82. 
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Fig. 5. — Initial bulge radius rbuige and initial bulge mass Mb u i ge as a function of circular velocity 
computed for the four halos shown in Figure 1 assuming z v i r { a \ = 3. 
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Fig. 6. — Disk and bulge infall rates vs. redshift for four DM halos with V^i rc = 200 km s and 
^viriai = 3 and r bu i ge = 4 kpc, A spin = 0.03. 
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Fig. 7. — Bulge blowout conditions for Type I and Type II supernova and a super-massive black 
hole. 



